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1. Introduction 

The understanding of the physics of cool white dwarfs (WD) bears im- 
portant consequences for Galactic evolution and cosmological implications. 
The observed cutoff in the disk WD luminosity function (WDLF) yields 
the determination of the age of the Galactic disk, as suggested initially 
by Winget et al. (1987). The recent microlensing observations toward the 
LMC (Alcock et al., 1996) suggest that WDs might provide a substantial 
fraction of the halo dark matter (Chabrier, Segretain & Mera, 1996; Adams 
&: Laughlin, 1996). The correct analysis of these applications implies a cor- 
rect WD cooling theory and reliable photometric predictions, which in turn 
require accurate interior and atmosphere models. Important improvement 
in this latter domain has been accomplished recently by Bergeron, Saumon 
& Wesemael (1995) and Bergeron, Wesemael & Beauchamp (1995), which 
yields the determination of photometric color indices and bolometric cor- 
rections down to 4000 K (see Leggett, these proceedings). In this paper, we 
review the most recent improvement in WD interior and cooling theory. 

2. Internal structure. Cooling theory 

The theory of WD cooling was first outlined by Mestel & Ruderman (1967) 
who identified the dichotomic properties of WDs, where degenerate elec- 
trons provide the pressure support but do not contribute significantly to the 
heat capacity (cv/ks = n'^ /2{T/Tp) for a fermion gas, where cy is the spe- 
cific heat per unit mass and Tp ~ 10^ K is the electron Fermi temperature 
for WD central densities) while the ions contribute negligible pressure but 
provide most of the thermal energy. Abrikosov (1960) and Salpeter (1961) 
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independently pointed out the possible onset of crystallization in cool WDs. 
VanHorn (1968) first developed a consistent theory of WD crystallization 
and Lamb & VanHorn (1975) first calculated the evolution of a pure carbon 
crystallizing WD. These calculations have been extended to C/0 mixtures 
by Wood (1992) who examined extensively the importance of the various 
parameters entering WD evolution, e.g. the core composition and the atmo- 
sphere structure. A further significant breakthrough in WD cooling theory 
is due to Stevenson (1980) who first pointed out the importance of the 
crystallization diagram in a two-component (e.g. C/0) plasma and of the 
difference of chemical composition in the fluid and in the solid phase. This 
motivated numerous calculations for the characterization of the phase dia- 
gram (Barrat, Hansen Sz Mochkovitch, 1988; Ichimaru, lyetomi Sz Ogata, 
1988), the effect on WD cooling (Mochkovitch, 1983) and on the WDLF 
(Garcia-Berro, Hcrnanz, Mochkovitch & Iscrn, 1988). More recently Seg- 
retain &: Chabrier (1993) characterized the evolution of the crystallization 
diagram of stellar plasmas for arbitrary binary mixtures as a function of the 
charge ratio. The chemical differentiation at crystallization calculated with 
these diagrams was shown to produce an extra source of energy AE in the 
WD (Segretain ct al., 1994), which in turn leads to a substantial increase 
in the age of crystallized WDs for a given luminosity {At oc AE/L), an 
important issue for a correct determination of the age of the Galactic disk 
(Hernanz al., 1994). 

In this section, we derive an analytical theory for the evolution of cool, 
crystallizing WDs, based on first principles of thermodynamics (Landau 
& Lifschitz, 1980), aimed at describing the main physical effects in terms 
of simple physics (see also Isern et al., 1997). These calculations include 
all the afore-mentioned processes and yield a reasonable estimate of the 
gravitational energy release and time delay induced by chemical fractiona- 
tion at crystallization, a question of strong debate among the community. 
This provides useful guidelines to verify the validity of complete, numerical 
calculations. 

The first laws of thermodynamics yield for white dwarf cooling : 



L = — —am — / ej/Cim = — / am — / e^am (1) 

Jo at Jo Jo at Jo 

where e^, is the neutrino rate and dq/dt is the heat rate per unit mass 
along the change of an equilibrium state. d,u and arc the change of 
specific internal energy and gravitational energy, respectively. The first one 
reads : 



dP 

du = cvdT + {{T—)y - P)dv + Auarys + [T,i{iiidNi)i + ^i{nidNi), 
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= Tds - Pdv + TAs + [^i{nidNi)i + J:i{piidNi)s]^^^ (2) 

The Hi denote the chemical potentials, s the specific entropy and v = l/p 
is the volume per unit mass. The dNi are the variations of carbon and oxy- 
gen nuclei in the fluid and the solid phase due to the change of composition 
at crystallization. They can be calculated with the thermodynamics lever 
rule. When the central crystal grows-*^, there is a thin C-enriched surround- 
ing fluid layer where locally the variation of nuclei d(6N^ ) is not equal to 
—dNi^ (no local mass conservation) ^. How this carbon excess in the fluid 
is redistributed homogeneously will be examined below. No variation of 
composition yields of course 5Ni = and thus no extra internal energy. 

The first two terms on the r.h.s. of (2) have been recognized originally 
by Mestel & Ruderman (1967), the third term was first introduced by 
Van Horn (1968) and is the crystallization latent heat I = —Aucrys = 

T{Ssol - Sliq). 

The condition of hydrostatic equilibrium yields for the variation of grav- 
itational energy : 



i-M p 

Sn = / PSvSm ^< — >M (3) 
Jo p 

since the electron pressure largely dominates the ionic pressure. 
Eqns. (l)-(3) can be rewritten : 



, , dT ^ f^,^dP, dv , ,dms , dms 

where du = Si(/ iJ,id{6Ni)) and nis is the mass crystallized. Note that 
the last two terms in eqns. (2) and (4) are evaluated at constant volume 
and thus do not stem from a contraction work, but from the change of 
composition at crystallization. 

The contribution of the first three terms of eqn. (4) can be estimated 
easily, as done initially by Mestel & Ruderman (1967) and Lamb &; Van 
Horn (1975) : 

^Notc that this fractionation process is described sometimes in the literature as drown- 
ing O-flakes. This is an erroneous picture, based on Stevenson's old eutectic diagram, 
which corresponds to an inhomogeneous sohd. The correct spindle diagram yields a ho- 
mogeneous solid. It is easy to verify that the solid is denser than the liquid, so that the 
crystal grows at the center of the star. Were WDs made of water, the story would be 
different because of the volume expansion at crystallization 

^This is similar to e.g. silicium deposition for semi-conductor devices in a sili- 
cium-l-impurity liquid, which yields a silicium concentration gradient and eventually a 
silicium-rich layer. We have in fact a distillation process in the WD. 
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I'M 

/ cydTdm^ from the virial theorem 
Jo 

fM fip rM fip. pth 

/ dvdm ^ / [T^ + oiT/TpfUvdm f=i< > M 

Jo dT Jo dT p 

(5) 

where P/'* = pTZT/jiQ (where jjlq is the mean ionic molecular weigth) is 
the thermal (non electrostatic) ionic pressure. 

Equations (3) and (5) yield SUgrav ~ > SO,. WD characteristic cen- 
tral density p ^ 10^ g.cm~^ and central temperature T rsj IQ^ K yield 
Pl^/Pe ~ 10~^/Z. Thus only a negligible fraction of the energy due to 
gravitational contraction is radiated. Most of the work is expended in rais- 
ing the electron Fermi energy, as first noted by Lamb & Van Horn (1975). 
The latent heat contribution can be estimated from the differences between 
the solid and liquid ionic entropies (VanHorn, 1968) : I ~ —kTc/AH, where 
Tc is the crystallization temperature (~ 3 x 10^ for C/0, see Segretain & 
Chabrier, 1993), A is the mean atomic mass and H = 1 a.m.u.= 1.66x10"^^ 
g. This yields an energy release Uiatentheat ~ 10^'^ erg ~ 10"^ $7, where 
n = GM^/R ~ 10^^ erg is the WD gravitational energy. The negative sign 
indicates that the energy is emitted at crystallization. 

The last term can be estimated as follows : 



SUth = 



= (^^v^tSX ^ {^)y,T^X ^ An, (6) 

where X is the mass fraction of one of the components (say carbon), 
AX = Xi — Xg and Aui = — Ui^ is the difference of Madelung energy 
in the C/0 plasma between the fluid and the solid phase^. The Madelung 
energy (per unit mass) of the mixture reads : 



5/3 5/3 

u,/kT = aSfe^Ffe = aTe{X^ + (1 - X)^} (7) 

-^k ^1 ^2 

where a denotes the Madelung constant in the fluid or in the solid phase 
{as = — 0.9, a; = —0.899), Fg = e^/a^kT (og is the mean inter-electronic 
distance) and the index k denotes each ionic species (C^"*", 0^+). This yields: 

^In fact 5u is only a fraction of Aw, since it stems from the difference w.r.t. to the av- 
erage energy over the C-enriched layer. This is not consequential for the present estimate 
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(8) 



J 6u drus < APi > Ms 



APj Pj Ms 
Pi Pe M 



(9) 



S^l <Pe> M 



Note that Pi is now the ionic electrostatic pressure. For Zi = 6, Z2 = 

8, xi=.T2=l/2, wc get Pi/Pe ~ 10"2 - iQ-i, APi/P, ~ 2Ax{zl^^ - 

Z2^)/{Zi^^ + ^2^^) 0.5 X Ax, where Aa; ~ 0.1 - 0.3 is the differ- 
ence of carbon number concentration between the solid and the fluid phase 
at crystallization (Segretain & Chabrier, 1993). With M^/M ~ 0.1, this 
yields: 



in agreement with the detailed numerical calculations (see Figure 5 of 
Segretain et al., 1994). 

Chemical differentiation at crystallization thus provides an additional 
source of energy wich remains much smaller than the gravitational energy. 
But, as shown below, the release of this quantity at a low-luminosity phase 
of the evolution has a significant effect upon the lifetime of the star at these 
stages. 



With the afore-mentioned values we get At ~ 5 x 10^ yr at the begining 
of crystallization, Ms/M ~ 10% and L = 10~^-^ Lq and At ~ 2 x 10*^ yr at 
L = W-^-^Lq, Ms/M ^ 80%, the observed cutoff luminosity. These simple 
calculations show the importance of the time delay induced by chemical 
fractionation at crystallization for cool (faint) WDs, even though the cor- 
responding energy is small compared to the binding energy. 

An other issue concerns the redistribution of the excess of carbon in 
the fluid at crystallization since the solid core is 0-enriched, as obtained 
from the phase diagram. Since the fluid C-enriched layer around the crys- 
tal is lighter than the surrounding medium, a Rayleigh- Taylor instability 
develops locally, due to the variation of molecular weight. This problem 
has been considered in detail by Mochkovitch (1983) who showed that the 



su/sn ~ 10"^ - 10"^ 



(10) 




(11) 
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typical crystallization time (for a 0.6Mq WD, the crystallization velocity at 
the begining of crystallization is Vc ~ 10~^ Mwd/7x 10^ yr ~ 10"^ crn.s"^) 
is significantly larger than the convection time so that the liquid is likely 
to be rehomogeneized rapidly as crystallization goes on. 

3. Galactic implications 

The first application of these calculations concerns the age of the faintest 
WD ever observed, 5*0 439 — 26 (Ruiz et al., 1995). The trigonometric- 
parallax determination of this object yields an absolute magnitude My = 
17 A ± 0.3, Mhoi = 17.1 ± 0.1 (Bergeron et al., 1997) about 1 mag faint- 
ward of the observed cut-off of the WDLF of the Galactic disk (Liebert et 
al., 1988). The location of this object in the My ys V — I diagram, and a 
comparison with the photometric sequences of cool white dwarfs recently 
derived by Bergeron et al. (1995) yields the interpretation that it is a cool 
{Teff ~ 4500 K), massive WD (m ~ 1.2 Mq). Comparison with (pure car- 
bon) evolutionary models of Wood (1992) yields an age determination for 
this WD, t ~ 6.5 Gyr, substantially below the lower limit for the age of 
the Galactic disk determined by detailed WD cooling theory (Hernanz et 
al., 1994). The Wood sequences do not include the afore- mentioned release 
of gravitational energy due to C/0 differentiation at crystallization. Fig- 
ure 1 displays several isochrones of a C/0 WD for several masses (WDs 
more massive than 1.2 Mq have a different 0/Ne/Mg internal composition), 
obtained with a WD cooling theory including (solid line) and neglecting 
(dashed line) the fractionation process. The dotted lines display the lumi- 
nosity of ESO 439-26 and its mass, assuming either a H-rich (~ 1.1 Mq) 
or a He-rich (~ 1.2 Mq) atmosphere (sec Ruiz et al., 1995). As shown, 
ESO 439-26 is compatible with an age t'^ 10 Gyr, in good agreement with 
the most recent determination of the age of the disk, whereas neglecting 
differentiation would yield t <9 Gyr. 

An other important galactic application is the age of the disk inferred 
from the comparison between the observed and the theoretical WDLF. As 
shown by Hernanz et al. (1994) and by the present analysis, chemical frac- 
tionation yields 1 to 2 Gyr older ages (i.e. 10-20% increase), depending on 
the initial C/0 profile, w.r.t. estimates which do not include this process. 
Preliminary calculations along the present lines based on the recently ob- 
served WDLF (Oswalt et al., 1996) yield an estimate for the age of the 
cutoff at log L/Lq = —4.5, i.e. an age for the Galactic disk, to ~ 12 Gyr. 

4. Uncertainties in the theory 

An important uncertainty in the models is the initial C/0 composition 
and stratification. If the initial WD is already stratified, with oxygen ac- 
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Figure 1. Isochrones for different masses of crystallized WDs (logL/Le < —3.5) in a 
mass- luminosity diagram. The solid lines include C/O differentiation at crystallization 
whereas the dashed lines do not. The isochrones are the same in both cases and go 
from t = 7 Gyr to t = 15 Gyr from left to right. The vertical dotted line is the inferred 
luminosity for £S'0 439 — 26 (Bergeron et al., 1997), while the horizontal dotted lines are 
the inferred mass for a H-rich (1.1 Mq) and He-rich (1.2 Mq) atmosphere. Crystallization 
starts earlier for massive (and thus rapidly evolving) WDs, which then enter the rapid 
Debye cooling regime, which causes the bending in the isochrones. 



cumulated near the center, as suggested by Mazzitelli &: D'Antona (1986), 
the energy release by differentiation, and the related time delay arc reduced 
(see Segretain et al, 1994, Tables 3-4). The other main source of uncertainty 
concerns the model atmospheres and most importantly for the age of the 
disk, the thickness of the helium layer wich regulates the emergent heat 
flux. The thicker the layer the younger the age (Wood, 1992). These uncer- 
tainties translate in a ~ 1 — 2 Gyr uncertainty on the age of the disk. Futher 
progress in this direction is certainly essential for a better determination of 
the cooling and the observational properties of cool WDs. 

5. Conclusion 

We have shown in this review that substantial improvement in the theory 
of cool WDs has been accomplished within the past few years for the atmo- 
sphere as well as for the internal structure. We have shown that chemical 
fractionation at crystallization, although it liberates a negligible amount of 
gravitational energy, modifies substantially the cooling history of the star 
and yields an important time-delay for faint WDs. This process cannot 
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be ignored in accurate WD cooling theory. As noted by Mochkovitch (pri- 
vate communication), although the importance of crystallization of alloys 
in stellar plasmas (WDs in particular) is still strongly debated in the astro- 
physical community, it has been recognized long ago in geophysics (see e.g. 
Loper, 1984; Buffett et al., 1992). Although the nature of the plasma is dif- 
ferent, the physics of the process (thermodynamics and energy transport) 
is exactly the same. Important uncertainties remain essentially in the exact 
determination of the initial C/0 profile in the star and in the structure of 
the hydrogen and/or helium outer envelope. 
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DISCUSSION 

ROBERT KURUCZ: What about the effect of rotation on solidification? 

GILLES CHABRIER: It is not clear how rotation will affect the 
crystal growth, but there is no reason for it to prevent 
crystallization, which must occur below a certain temperature. The 
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effect of rotation on the C-redistribution in the surrounding layer 
has been considered in detail by Mochkovitch (1983) . 



